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ABSTRACT 

We explore the development of structures in molecular gas in the Milky Way by applying the 
analysis of the brightness distribution function (BDF) and the brightness distribution index (BDI) 
in the archival data from the Boston University-Five College Radio Astronomy Observatory ^"^CO 
J = 1-0 Galactic Ring Survey. The BDI measures the fractional contribution of spatially confined 
bright molecular emission over faint emission extended over large areas. This relative quantity is 
largely independent of the amount of molecular gas and of any conventional, pre-conceived structures, 
such as cores, clumps, or giant molecular clouds. The structured molecular gas traced by higher BDI 
is located continuously along the spiral arms in the Milky Way in the longitude-velocity diagram. 
This clearly indicates that molecular gas changes its structure as it flows through the spiral arms. 
Although the high-BDI gas generally coincides with H II regions, there is also some high-BDI gas with 
no/little signature of ongoing star formation. These results support a possible evolutionary sequence 
in which unstructured, diffuse gas transforms itself into a structured state on encountering the spiral 
arms, followed by star formation and an eventual return to the unstructured state after the spiral arm 
passage. 

Subject headings: Galaxy: disk — ISM: molecules — radio lines: ISM 



1. INTRODUCTION 

A stellar spiral potential causes a large amount of gas 
to ac c umulate along the spiral arms (e.g., Wada & Koda 



200l iDobbs et all 12001 iShettv fc 6strikei1 l2006t iWada 
20081) . The concentration and compression of the molec- 
ular gas in spiral shocks leads to active star form ation 
(;e.g.. lRoberta[T96llRaIdl[T99llKnapen et al.|[T996t ). and 
indeed, most H II regions are found along the narrow dust 
lane s and molecular spiral arms in external spiral galax- 
ies (iRand et al.l 1 199a iScoville et all 120011 : lEgusa et al.) 
120111) ^ The physical conditions of the molecular gas 
also appear to change upon the entry into a spiral arm. 
Analyses of the CO line ratios suggest that the molec- 
ular gas becomes denser and/or warmer in spiral arms 
(ISakamot o et all Il997t iSempere fc Garci'a-Burilld 119971: 
pTosaki et al.l 120021 ^ A remarkable contrast of the CO 
J = 2-1/ J = 1-0 ratio is detected between the spiral 
arms and inter-arm regions in the grand-design spiral 
galaxy M51, suggesting an evolution of the physical con- 
ditions of the gas from the inter-arm regions, through 
the spiral arras, an d back into the next inter-arm regions 
()Koda et al.l 120121 ). These recent developments, how- 
ever, are based on low-resolution observations of external 
galaxies (often > 1 kpc) . Questions remain as to how the 
structures of molec ular gas, especially o n the scale of star 
formation (~ 1 pc; iLada fc Ladal 120031 ) . evolve through 
spiral arm passages, and how the evolution of the gas 
structures results in star formation. 

The Milky Way is the nearest spiral galaxy in which 
we can resolve pc-sized molecular gas structures (i.e., the 



scale directly relevant to sta r and cluster formation) . Re- 
cently, ,Sawada_et.M] (HnH) (hereafter iPaper I) have pre- 
sented high-resolution ^^CO J = 1-0 and ^^CO J = 1-0 
images of a 0?8 x 0?8 area through the disk of the Milky 
Way at I = 38°, and have shown that the spatial struc- 
ture of the molecular gas varies distinctly between the 
arm and inter- arm regions. They found that the bright 
and spatially confined emission (typical scale of several 
pc) is predominantly seen in the Sgr arm, while the 
fainter, extended emission dominates in the inter-arm re- 
gions^ 
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iPaper II quantified this structural evolution, by in- 
troducing the Brightness Distribution Function (BDF) 
and Brightness Distribution Index (BDI). In particu- 
lar, the BDI is a measure of the fractional contribu- 
tion of spatially confined bright molecular emission over 
faint emission extended over large areas, and can be de- 
rived directly from velocity channel maps without pre- 
assumption of structures. Their analysis, however, is 
limited to this relatively small region, and it is impor- 
tant to determine the BDF and BDI over the Milky Way 
disk to reveal whether this important structural evolu- 
tion is a general, fundamental process that leads to star 
formation across spiral arms. 

In th is Letter, we apply the analysis developed in 
IPaper II to archival "CO J = 1-0 line emission data 
covering the majority of the Galactic plane in the first 
Galactic quadrant. 

2. DATA AND ANALYSIS 

2.1. GRS Data 

We use the ^■^CO J — 1-0 data from the Boston 
University-Five College Ra dio Astronomy Obse rvatory 
Galactic Ring Survey (GRS. [Jackson et"alll2006D . which 
covers a large area of the Milky Way disk, 18° < / < 56°, 
~1° ^ ^ ^ -1-1°. The data are converted from the an- 
tenna temperature into the main-beam temperature 
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Tmb, by dividing by the main-beam efficiency, 0.48. 
The spatial resolution and sampling are 46" and 22", re- 
spectively. Although the velocity sampling of the original 
data is 0.21 km s~^, we take a running mean with a width 
of 2 km s^^ sampled every 1 km s~^, in order to improve 
the signal-to-noise ratio. The resultant velocity resolu- 
tion and sampling are similar to those used in iPaper II 
(2.7 and 1.3 km s~^, respectively). The BDI is calculated 
in regions of Al x Ab = 1° x 2° and Av = 1 km s^^. 

2.2. The Brightness Distribution Index 

Remarkable variations in gas structure are visually ap- 
parent throughout velocity channel maps. Figure [ija) 
andlllb) show the maps of the two characteristic regions, 
H-1 and L-1 (their selections and properties are discussed 
in Section [22]). In H-1, the emission shows high contrast 
in the spatial structure: i.e., there exist many compact 
and bright (> 4 K) structures, as well as diffuse and faint 
(^ 1 K) emission. In L-1, on the other hand, the con- 
trast is much lower than that in H-1: the diffuse and 
faint emission dominates, while there are only very few 
compact, bright structures. 

The BDF is defined in a certain volume in l-b-v space 
as the fraction of the "pixels" with brightness between 
T and T + dT (i.e., the histogram of the brightness). It 
illustrates the differences in the spatial structure of the 
emission noted above. Figure [Ijc) displays the BDFs in 
the regions H-1 and L-1. The presence of the compact 
and bright structures in H-1 appears as a long tail of the 
BDF toward high brightness, while the BDF in L-1 drops 
st eeply. 

'Paper I introduced the BDI as a way to characterize 
the BDF in a single number, so that it can be correlated 
with other parameters. It is defined as the flux ratio of 
the bright emission to faint emission: 



BDI = log 



7jr.B(T)dT\ 
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where B{T) denotes the BDF; Tq, Ti, T2, and T3 are the 
brightness thresholds; and T[i] is the brightness of the 
ith pixel in the l-b-v space. I n this Le tter we adopt the 
same brig htness thresholds as IPaper J . (Tq, Ti, Ts, Tg) = 
(1, 1.5, 4, cxo) [K]. A high BDI denotes dominance of com- 
pact, bright structures. The BDIs in the regions H-1 and 
L-1 are —0.63 and —2.15, respectively. 

The BDI is, by its definition, (1) free from any assump- 
tions about structures in the gas, such as giant molecu- 
lar clouds (GMCs) and cores/clumps within GMCs; and 
(2) independent of the amount of gas (i.e., an intensive 
property). The BDI is meaningful when it is defined 
with a number of pixels sufficiently large to represent the 
statistical characteristics of the volume of interest, and 
with a spatial resolution better than the typical scale of 
bright, compact structures (several pc: Paper I) so that 
such structures are not smeared out. One degree corre- 
sponds to « 50, 90, and 140 pc at distances of 3, 5, and 
8 kpc, respectively. Therefore, each area in the l-b plane 
(1° X 2°) that we analyze here is larger than the ^tradi- 
tional GMC whose average size is ~ 40 pc (jScoville et al.l 



Il987f ). The spatial resolution of 46" corresponds to 0.7, 
1.1, and 1.8 pc at 3, 5, and 8 kpc, respectively. 

3. RESULTS 

Figure [2l^a) and (b) present the longitude-velocity {l- 
v) distribution of ^■^CO intensity and the BDI, respec- 
tively. Figure [21(c) displays the line intensity and the 
BDI simultaneously: i.e., the brightness represents the 
intensity [same as (a)] and is colored by the BDI value 
[same as (b)]. The H II region s which are located inside 
the sky coverage of the GRS ([Anderson fc Banial[200a 
are overlaid on the l-v diagram of the BDI in Figure 
[21(b). The white polygons on Figure [2jb) are the regions 
which are discussed in Section [321 



3.1. Enhanced BDI along Spiral Arms 
There is a band of high BDI (> —1) from {l,v) « 
(20°, 30 kms^i) to (50°, 60 km s^^) in the l-v dia- 
gram. This high-BDI band coincides with a band of 
H II regions which was found in space and velocity us- 
ing radio recombination lines and has been commonly 
identified as a spiral arm of the Milky Way, the Sgr 
arm (e.g.. iMezged 119701: IGeorgelin fc GeorgelinI 119761 : 



iLockmanl 119791: iDownes et all 1 198 J). Chains of H II 
regions are inde ed found along spira l arms in external 
galaxies as well ([Scoville et al.|[200l[ ). It is noteworthy 
that the region ss (30°, 40 km s""'^) on the band shows 
high BDI, even though the ^'^CO intensity is low. The 
structured molecular gas traced by the higher BDI re- 
veals the continuation of the spiral arm ind ependent of 
the total intensity (amount) of the gas. In [Paper II we 
found that the BDI is high at the spiral arm velocities in 
a 0?8 X 0?8 field at Z « 38°. The present result supports 
our finding with analysis of data which covers a much 
wider region. The Set arm also exists in the coverage 
of the GRS, from {l,v) « (20°, 60 kms-^) to the Set 
tangent at (30°, 100 km s~^), and shows high BDI. 

The high-BDI gas distributed along the loci of the spi- 
ral arms corresponds to the "hot" or "warm" molecular 
clouds (i.e., bright in the ^^CO line) which also show a 
confined distri bution in the l-v domain compared with 
fainter clouds (ISanders et al.lll985l : [Solomon et al.l [19851 : 
IScoville et al.lll987[) . 

3.2. Individual Regions 

The polygons on Figure [^Jb) are example of regions 
with three ranges of BDI: low (L-1 and L-2; < —2), mod- 
erately high (M-1 and M-2; ~ -1), and high (H-1 and 
H-2; > —0.5), as well as an inter-arm region (I-l). These 
are chosen as representative regions, based on inspection 
of the full data cube, and with environments well char- 
acterized in the literature. 

3.2.1. Massive Complexes with Low BDI 

We find two complexes which have low BDI in spite 
of high intensity. One of them, L-1, is likely to be 
located at the near side of I ~ 40° (distance = 2.2 
kpc. IDame fc ThaddeusI [l985lfl . and was shown in Fig- 
ure [Hb). The other (L-2) is at the tangent point of 



Although IDame & Thaddeu! 
as a W50 cloud and IHuang'er^ 



J1985I) called this complex 
l |198 3l) discussed the physi- 
cal association betwee n the cloud and W50 SNR, a recent study 
HLockman et al.ll2007l 1 suggested that W50 is located in the far 
background. 
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r„,(''coy=i-o) [K] 

Fig. 1. — The velocity channel maps of the regions (a) H-1 and (b) L-1. The averaged Tmb in a 5 km s~^ channel (50-55 and 30-35 
km , respectively) is shown in the same linear scale, (c) The BDFs in these regions. 
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Fig. 2. — (a) The longitude— velocity distribution of the ^"^CO intensity averaged over Al X Ab = 1° X 2°. The loci of the Sgr arm 
(IS andcrs ct al. 1985) and the Set arm (Dame & Thaddcus 2011) are overlaid, (b) The BDI (pseudo color), the H II regions (white circles, 
ITndcrson & Bania 2009), and the regions to be discussed (white polygons). The region where the BDI cannot be defined (i.e., no emission 
brighter than To) and/or the ^^CO intensity is less than 0.01 K is represented by gray, (c) Images (a) and (b) combined: i.e., the brightness 
and color represent the ^'^CO intensity and the BDI, respectively. 
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Fig. 3.— The BDFs in the regions I-l (black), L-1 (blue), L-2 
(cyan), M-1 (green), M-2 (light green), H-1 (red), and H-2 (pink). 

/ ~ 38° (6.7 kpc lPaper it ). They are located outside the 
l-v loci where the Set and Sgr arms exist, and therefore, 
are supposedly in the inter-arm regions. The BDFs in 
these regions (Figure [Sj drops steeply toward high inten- 
sity, whi ch is similar to that we found in a part of L-2 
ijPaper J ). We also derived the BDF in the inter-arm re- 
gion LI whose intensity is much lower than those of L-1 
and L-2. The slope of the BDF in LI is close to those in 
L-1 and L-2 at Tmb > 1 K. The offset of the BDFs (i.e., 
LI is lower by c± 1) is attributed to the lower filling factor 
in l-b-v volume in LI by about an order of magnitude. 

L-1 and L-2 are distinguished from other inter- 
arm regions by their high intensity (i.e., concentra- 
tion of molecular content). Assuming local thermo- 
dynamic equilibrium at an excitation temperature of 
10 K and the abundance ratios of [^^C0]/[H2] = 



10-4 and fi^COl/r^COl = 50 ( 


Langer & Penziad 119931: 


iSavaee et al.ll2002l iMilam et al.l 


2005|). we estimate their 



masses to be 0.26 and 4.4 x IO^Mq, respectively. 
However, a substantial part of L-1 is extended be- 
yond the latitude coverage of the GRS, |6| < 1°. 
iDame fc"^h addeus (1985) estimated the mass to be 
2x 10^ Mp, base d on th eir wide-field ^^CO J = 1-0 obser- 
vations. iShand (|1972D derived the atomic mass density 
(~ 3 X IO^Mq kpc"^) and the azimuthal extent (~ 20°) 
of the H I feature corresponding to L-2. These values 
yield the atomic mass of ~ 6 x 10^ Mq, which is compa- 
rable to the molecu lar mass. 

The catalog of lAnderson fc Banial ()2009l ) lists no 
H II regions in L-1 , and only a few in L-2. 
iMaddalena fc ThaddeusI ()1985[ ) noted that the complex 
L-1 might be an object similar to the molecular cloud 
they found (Maddalena's cloud) due to the lack of signs 
of star formation despite their large masses. It is interest- 
ing that such massive clouds (~ 10^ Mq, w hich is near the 
uppe r end of the mass spectrum of GMCs: [Solomon et al.l 
|1987[) have low BDI (lack structure) and show little ev- 
idence of star formation. The regio n L-2 has been pro - 
posed to be an inter-arm spur (e.g., IDame et al.l 119861 ) . 
This might be a counterpart of the spurs seen in exter- 
nal g alaxies (e.g.. lAdler et al.|[l992HGarcia-Burillo et al.l 
\mm . which are suggested to be remnants of massive 
GMCs a s they are shredde d and expand after spiral arm 
passage (|Koda et al.ll2009D . 



3.2.2. High-BDI (structured) Gas and Star Formation 

The l-v distribution of H II regions in the first Galactic 
quadrant follow s narrow patter ns which trace the Set and 
Sgr arms (e.g.. iLockmanl [19791) . On these patterns, the 
BDI tends to be high (> — 1), as men tioned in Section 
13.11 It supports our finding in [Paper II that high BDI is 
seen at the velocities where H II regions exist. Hence, in 
a global point of view, high-BDI (structured) molecular 
gas is coincident with star-forming regions. Two inter- 
pretations are possible for this spatial coincidence: the 
high-BDI molecular gas could be either a result of star 
formation, or the cause. 

There are indeed some moderately-high-BDI regions 
associated with few H II regions, which may indicate 
that the development of structured gas precedes star 
formation. The regions M-1 and M-2 on the Set arm 
show moderately high (~ —1) BDI and high ^■^CO in- 
tensity, although the density of H II regions within them 
is low compared w ith the other parts of the Set arm. 
IDame et al.l (|1986l ) also pointed out that M-2 is an un- 
usual complex in that it lacks associated H II reg ions 
(only one in the catalog of lAnderson fc Baniall2009l) . 

4. DISCUSSION 

The above results may indicate the following scenario 
for the evolution of the gas through the spiral arms. Faint 
and diffuse (not structured) molecular gas, whose proto- 
types are the I-l, L-1, and L-2 regions, exist in the inter- 
arm regions. Entering the spiral arms, the gas develops 
bright and compact structures which become pre-star- 
forming complexes like M-1 and M-2. When stars form 
from these structures, the BDI becomes very high. The 
regions H-1 and H-2 may be in this stage. Then, once 
the gas leaves the arms, it returns to a diffuse state. 

Comparison of the BDFs in individual regions gives 
us a clue to understanding how star formation affects 
BDF/BDI. The regions H-1 (W44) and H-2 (W43) are 
prototypes of active star-forming regions (i.e., associ- 
ated with a number of H II regions) and show very high 
(> -0.5) BDI. The BDFs in these regions (Figure [2]) are 
characterized by long tail toward the high brightness (> 8 
K), which distinguish them from the moderately-high- 
BDI regions M-1 and M-2 (presumably pre-star- forming 
complexes). While M-1 and M-2 show a similar level of 
emission at Tmb ~ 4 K as H-1 and H-2, they lack emis- 
sion at Tmb ^ 8 K. This suggests that star- forming ac- 
tivity produces an excess in the BDF at very high bright- 
ness. 

From the comparison of the BDFs in ^^CO and ^^CO in 
[Paper II (Figures 7 and 8 therein), the emission brighter 
than 8 K in ^^CQ J = 1-0 corresponds to Tmb > 20 K in 
^^CO J = 1-0, which indicates the kinetic temperature 
Tk > 25 K if the ^^CO line is thermalized and opti- 
cally thick. This temperature is higher than the kinetic 
temperature of the bulk molecular gas in GMCs, 10-20 
K, corresponding to Tmb ~ 3-6 K in ^^CO J = 1-0. 
Therefore the long tail of the BDF at Tmb > 8 K in the 
high-BDI regions is likely the consequence of the intense 
heating around star- forming regions, e.g., in photodisso- 
ciation regions (iTielens fc Hollenbach|[l985[) . 

In order to separate the two possibilities, i.e., the high- 
BDI gas as result or cause of star formation, we exclude 
the gas under direct influence of massive star formation. 
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Fig. 4. — (a) The ^'^CO intensity (brightness) and BDI (color), the same as Figure[2{c). (b) Same as (a), but the radii of < 10 pc from 
ultra-compact and compact H II regions and < 25 pc from diffuse H II regions are excluded at the velocity channels within < 10 km s~^ 
from the H II regions. For the H II regions whose distances are not determined, the near kinematic distances are assumed, (c) The BDF 
in the region H-1, before and after excluding the gas around H II regions. 

along the bands of H II regions (i.e., the Sgr and Set 
arms). Molecular gas in the spiral arms is struc- 
tured, hosting relatively abundant compact con- 
centrations. 



Wc mask the radii of < 10 pc from the ultra-compact 
and compact H II regions an d < 25 pc from the diffus e 
H II regions in the catalog of I Anderson fc Banial ()2009[ ). 
Figure ID^b) displays the l-v diagram after the mask. 
Compared with the original BDI (Figure Hl^a)), high-BDI 
gas mostly turns into moderately-high-BDI gas, whereas 
moderately-high-BDI gas persists in the spiral arms. The 
BDF in the region H-1 before and after excluding the 
H II regions is shown in Figure lUJc). The BDF after 
excluding the H II regions is truncated at Tmb ~ 8 K, 
resembling the BDFs in the moderately-high-BDI regions 
M-1 and M-2. These results reinforce our argument that 
moderately-high-BDI gas formation precedes star forma- 
tion, while the high-TMB tail of the BDF in high-BDI 
regions is the consequence of star formation. This anal- 
ysis alone cannot exclude other possibilities that could 
enhance the BDI, such as due to potential uncatalogued 
H II regions and supernova remnants. More extensive 
studies are necessary to assess the impacts of their ef- 
fects. 

By applying the BDF/BDI analysis to the high spatial 
resolution data, the present study enables us to distin- 
guish between structured molecular gas with and with- 
out star formation, both of which were cl assified as the 
same pop ulation in low-resolution studies (I Sanders et al.l 
[1985; Sol omon et all 119851: IScoville et al.l 119871 see Se~ 
tion l3.ip . More studies of the BDF/BDI in nearby molec- 
ular clouds would provide us with better understanding 
of the relationship between the star-forming activity and 
the BDF/BDI. 

5. CONCLUSIONS 

In this Letter, we examined the structural variation 
of molecular gas in the first Galactic quadrant. We ap- 
pl ied the a nalyses of the BDF and the BDI developed 
in IPaper II to the ^^CO J = 1- emission line data from 
the GRS (| Jackson et al.l 120061) . The major findings are 
summarized as follows: 

1. The high BDI bands in the l-v diagram are found 



2. In the inter-arm regions, we found two examples of 
massive molecular gas complexes (~ 10^ Mq) which 
show low BDI and no/little signature of star for- 
mation despite their large masses. 

3. Molecular gas with high BDI usually has H II 
regions associated with it in space and velocity. 
There is also an interesting component of the 
molecular gas with moderately-high BDI which 
shows no/little signature of ongoing massive star 
formation (H II regions). 

4. These results suggest the structural evolution of 
molecular gas through spiral arms passages. The 
faint and mostly extended gas in the inter-arm re- 
gions develops bright and compact structures upon 
entering the spiral arms. Stars form in the compact 
structures of the gas. The gas then becomes diffuse 
as it leaves the arms. 

As we see in this Letter, the BDF and BDI of molecu- 
lar lines are new, useful tools to diagnose the properties 
of molecular gas, in terms of the spatial structures in the 
gas. The physics that underlie the observed evolution is 
still unknown, though it is likely related to the gas dy- 
namics around spiral arms. With the full commissioning 
of Atacama Large Millimeter/submillimeter Array, high- 
fidelity molecular line images of nearby galaxies with the 
critical several-pc spatial resolution is becoming feasible. 
Such images would enable us to make BDI maps for ex- 
ternal galaxies, which should open a door to a synthetic 
picture of the evolution of interstellar gas and star for- 
mation in various galactic environments. 

This publication makes use of molecular line data from 
the Boston University-FCRAO Galactic Ring Survey 
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